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Abstract 

In combustion environments, volatilization of Si0 2 to Si-O-H(g) species is 
a critical issue. Available thermochemical data for Si-O-H(g) species were used 
to calculate boundary layer controlled fluxes from SiO z . Calculated fluxes were 
compared to volatilization rates of Si0 2 scales grown on SiC which were 
measured in Part I of this paper. Calculated volatilization rates were also 
compared to those measured in synthetic combustion gas furnace tests. 
Probable vapor species were identified in both fuel-lean and fuel-rich combustion 
environments based on the observed pressure, temperature and velocity 
dependencies as well as the magnitude of the volatility rate. Water vapor is 
responsible for the degradation of Si0 2 in the fuel-lean environment. Silica 
volatility in fuel-lean combustion environments is attributed primarily to the 
formation of Si(OH) 4 (g) with a small contribution of SiO(OH) 2 (g). Reducing gases 
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such as H 2 and/or CO in combination with water vapor contribute to the 
degradation of Si0 2 in the fuel-rich environment. The model to describe silica 
volatility in a fuel-rich combustion environment gave a less satisfactory fit to the 
observed results. Nevertheless, it was concluded, given the known 
thermochemical data, that silica volatility in a fuel-rich combustion environment is 
best described by the formation of SiO(g) at one atmosphere total pressure and 
the formation of Si(OH W g). SiO(OH) 2 (g), and SiO(OH)(g) at higher pressures. 
Other Si-O-H(g) species, such as Si 2 (OH) 6 , may contribute to the volatility of 

Si0 2 under fuel-rich conditions, however, complete thermochemical data are 
unavailable at this time. 

Introduction 

SiC and other Si-based ceramic composites have been proposed for high 
temperature components in aircraft turbine engines. In Part I of this paper [1], 
high pressure burner rig (HPBR) testing of CVQ SiC in both fuel-lean and fuel- 
rich environments resulted in linear weight loss and surface recession rates as a 
result of silica volatility. The objective of this paper is to develop a chemical 
model for the volatilization of Si0 2 in complex combustion environments which 
describes the weight loss and surface recession observed for SiC in the HPBR, 
as well as other furnace tests of SiC in synthetic combustion environments. 

The combustion of Jet A Fuel, CH, . 9 , 85 , in air breathing engines gives rise 
to an environment which always contains N 2 , H 2 0, and C0 2 . In fuel-lean 
environments 0 2 is also found. On the other hand, fuel-rich environments 
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contain H 2 and CO but little O 2 . The gas composition of the combustion 
environment has been discussed in greater detail [2] and is shown in Figure 1 as 
a function of equivalence ratio, <J>. Equivalence ratio is a fuel-to-air ratio with total 
hydrocarbon content normalized to the amount of oxygen. At stoichiometry, <(>=1, 
combustion results in complete consumption of fuel and air. Regardless of fuel- 
to-air ratio, about 10% of the combustion gas is composed of water vapor. It will 
be demonstrated in this paper that water vapor is primarily responsible for the 
degradation of SiC in the combustion environment. 

Theory and Literature Review 

SiC is thermodynamically unstable in an oxidizing environment and forms 
an outer scale of S 1 O 2 . Because the Si02 forms a protective layer which grows 
at a slow rate, SiC has been proposed for use in high temperature oxidizing 
conditions such as combustion environments. In combustion environments 
containing O 2 , CO 2 , and H 2 0, SiC can oxidize by any or all of the following 
reactions: 

SiC + 3/2 02 (g) = Si02 + CO(g) (■) j 

SiC + 3C0 2 (g) = Si02 + 4CO(g) (2) 

SiC + 3H20(g) = Si02 + 3H 2 (g) + CO(g) ( 3 ) 

Water vapor has been found to be the primary oxidant based on a comparison of 
the oxidation rates of SiC in each gas [3]. 
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It has also been shown that in mixed oxidizing/reducing gases [4-6] such 
as H 2 0/H 2 or C0 2 /C0 mixtures, the silica scale can, in turn, be reduced by one 
of the following reactions to form volatile SiO(g): 

Si0 2 + H 2 (g) = SiO(g) + H 2 0(g) (4) 

Si0 2 + CO(g) = SiO(g) + C0 2 (g) (5) 

Similarly, in water vapor containing environments the silica scale may react to 
form volatile hydroxides or oxyhydroxides by one of the following reactions: 

Si0 2 + H 2 0(g) = SiO(OH) 2 (g) (6) 

Si0 2 + 2H 2 0(g) = Si(OH)4(g) (7) 

2Si0 2 + 3H 2 0(g) = Si 2 0(0H) 6 (g) (8) 

Si0 2 + V2 H 2 0(g) = SiO(OH)(g) + 1/4 0 2 (g) (9) 

2Si0 2 + 3H 2 0(g) = Si 2 (OH) 6 (g) -t- 1/2 0 2 (g) (10) 

Reaction 7 has been observed experimentally by a transpiration method [7] as 
well as atmospheric sampling mass spectrometry [8]. Reactions 6 and 9 have 
been observed by Knudsen effusion mass spectrometry [9,10]. Reaction 8 is 
thought to be important at pressures much higher than those found in turbine 


engines [11]. Reaction 10 has been described by Krikorian [11], but has not 
previously been considered important in water vapor containing environments. 
No experimental data are available for this molecule and estimated data are 
incomplete. 

In conditions, such as combustion environments, where both SiC 


oxidation and Si0 2 volatilization occur, paralinear kinetics are observed. The 
oxidation/volatilization kinetics have been modeled by Tedmon [12] for chromia 
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forming materials and for SiC by Opila & Hann [13]. Paralinear kinetics, 
expressed in terms of specific weight change, typical for SiC exposed in 50% 
H 2 O/50% 0 2 at 1200°C, are shown in Figure 2a. The overall weight change is 
due to the sum of the weight gain due to the growth of the scale and the weight 
loss due to volatilization of the silica. At long times, oxide growth occurs at the 
same rate that oxide volatilization occurs so that a constant oxide thickness is 
formed while a linear weight loss is observed. Figure 2b shows the 
corresponding dimensional change of SiC during paralinear oxidation. The rate 
constants used to plot these curves are given in terms of oxide thickness rather 
than weight change, but are equivalent to those used in Figure 2a, so that direct 
comparisons may be made. Note that even after a constant oxide thickness is 
achieved, SiC recession continues at a linear rate as silica volatilizes. Under 
conditions where the volatility rate is much greater than the oxidation rate, nearly 
linear weight loss and recession will be observed even at short times. The rate 
of SiC recession is thus controlled by the volatility rate of silica rather than the 
oxidation rate of SiC. 

The flux of volatile silicon species is limited by diffusion through a gaseous 
boundary layer. Boundary layer limited fluxes, J, can be calculated for a flat 
plate using the following equation[14]: 

J = 0.664 Re 1/2 Sc V3 ^ (11 ) 

where J is a mass flux, Re is the Reynolds number, Sc is the Schmidt number, D 
is the interdiffusion coefficient of the volatile species in the boundary layer 
combustion gas, p v is the concentration of the volatile species at the solid/gas 
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interface and L is the characteristic length of the test specimen parallel to the 

direction of the gas flow. Expanding the Reynolds and Schmidt numbers results 
in: 


J = 



(12) 


Here, p is the concentration of the boundary layer combustion gas, v is the linear 
gas velocity, and rj is the gas viscosity. Gas concentrations are calculated from 
the ideal gas law. The interdiffusion coefficient is calculated using the Chapman- 
Enskog Equation [15] using tabulated gas constant values [16,17]. Values for 
silicon hydroxide and oxyhydroxide species are approximated with values for 
corresponding fluoride molecules or similar oxide molecules. 

At higher gas velocities, laminar flow is expected to give way to turbulent 
conditions. In a smooth-walled tube of circular cross-section this transition 
occurs for Reynolds numbers, based on tube diameter Rep, greater than 2100 
[18]. A flat plate placed in turbulent flow will build up a laminar boundary layer 
near the leading edge. Flow along the plate will become turbulent at some 
distance downstream. This transition occurs at a Reynolds number, based on 
the distance from the leading edge Re x , of 2x1 0 5 to 5x10 s [17]. Under these 
conditions, the boundary layer of combustion gases is disturbed, allowing faster 
volatilization of the surface layers. In the limiting case, the boundary layer would 
be completely swept away and the volatilization rate would depend only on the 
vapor pressure of the volatile species, equivalent to volatilization into a vacuum. 
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This maximum possible mass flux can be calculated using the Langmuir 
Equation: 


, PM 1/Z 

(2rcRT) 1/2 ( 13) 

Here, the flux , J, is in units of mass per unit area per time, M is the molecular 
weight of the vapor species, R is the gas constant and T is the absolute 
temperature. Boundary layer calculations based on laminar flow, therefore, give 
minimum volatilization rates. 

In the cases discussed herein, boundary layer control dominates. 
Equations 11 and 12 can be recast in terms of known combustion parameters, 
such as gas velocity, total pressure, and partial pressures of the volatile species: 

v « 

^ ° C pl/2 ” P volatile (14) 

The temperature dependence of these fluxes is contained within the partial 
pressure of the volatile species. The partial pressures of the volatile species can 
in turn be written in terms of the combustion gas phase constituents: O z , H 2 0, 

H 2 , C0 2 , CO, and N 2 . Equilibrium constants can be written for reactions 4 
through 9. For example, 

[P(SI<OH)„)1 

[P(HaO)] 2 <15) 

From this simple expression it can be seen that the P(Si(OH) 4 ) depends on the 

square of the water vapor pressure. It follows, then, that for volatility due to 
Si(OH) 4 formation: 
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Jsi(OH) 4 ~ j^T-PHO) 2 (16) 

Mass fluxes obtained at a given temperature can thus be easily extrapolated to 
those at other gas pressures and velocities once the correct volatile formation 
reactions are identified. For each volatile species, the velocity exponent is 
predicted to be 1/2. However, the pressure and temperature dependence of 
each species will vary, and these dependencies are summarized in Table 1. 

The probable identity of the volatile species present in various combustion 
environments is made by comparing the overall magnitude of the fluxes, the 
pressure dependence, and the temperature dependence of the measured data 
to calculated values under the same conditions. Thermodynamic data can be 
used to predict vapor pressures of volatile species, which can in turn be used in 
Equations 11 and 12 to calculate boundary layer controlled fluxes of volatile 
species. While thermodynamic data for SiO(g) are readily available [19], data 
for hydroxides and oxyhydroxides are more limited. Hashimoto [7] has 
experimentally determined limited thermodynamic quantities for Si(OH) 4 (g) while 
Hildenbrand and Lau [9,10] have experimentally determined limited 
thermodynamic data for the SiO(OH)(g) and SiO(OH) 2 (g) species. It was, 
therefore, necessary to make use of several sets of estimated thermodynamic 
quantities [1 1 ,20] to calculate fluxes of the silicon oxyhydroxide species. The 
thermochemical data of Krikorian [11] are based on equilibrium constants for 
volatilization reactions from a compilation of studies as well as entropies 
calculated from spectral data and partition functions. The thermochemical data 



of Allendorf et al. [20] and Darling and Schlegel [21] are derived from ab initio 
calculations using Hartree-Fock theory. 

The aim of this paper is to use the available thermodynamic data and flux 
calculations to develop a self-consistent model for the volatilization of silica 
during paralinear oxidation of SiC in combustion environments. The 
experimental process [1,13,22] should be described with as accurate a chemical 
model as possible. Verification of this model will thereby enable the prediction 
of SiC recession under a variety of combustion conditions. 

PROCEDURE 

Experimental data from three sources were examined in light of possible 
Si02 volatilization/SiC recession reactions. These studies are the HPBR fuel-rich 
and fuel-lean conditions from Part I of this paper, synthetic fuel-lean furnace 
tests [13], and synthetic fuel-rich furnace tests [22]. Synthetic gas mixtures are 
premixed to reflect the components found in combustion gas environments. For 
example, the synthetic fuel lean mixture models the oxygen and water vapor 
found at <|>=0.6 and a total pressure of 5 atm, while the synthetic fuel-rich mixture 
models the combustion gas products at <^=1.5. The experimental parameters 
used in the calculations are shown in Table 2. 

The available thermodynamic data for volatile Si-O-H species have been 
adapted for use in ChemSage [23], a free energy minimization program. The 
free energy expressions for the species added to the ChemSage data base are 
shown in the Appendix. Data for SiO are from JANAF [19], data for Si(OH) 4 
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have been experimentally determined by Hashimoto [7], and two sets of 
estimated data for SiO(OH) and SiO(OH) 2 [11,20] were also used. Additional 
data for SiO(OH) by Darling and Schlegel [21] are referred to, but not used in 
the calculations. Si0 2 (cristobalite) plus the combustion gas composition at the 
experimental temperature and pressure were input to ChemSage [23], The 
calculations were each repeated with the two data sets specified in Table 3. 
Vapor pressures of SiO, Si(OH) 4 , SiO(OH), and SiO(OH) 2 were calculated in 
each case. The vapor pressure of Si 2 0(0H) 6 was found to be negligible under 
the conditions of this study. 

Boundary layer controlled fluxes of the volatile species were then 
calculated using the ChemSage vapor pressures in Equation 12. The boundary 
layer thicknesses, 5, for the low gas velocity experiments were calculated to vary 
between 0.69 and 0.75 cm for the synthetic fuel-lean furnace experiments [13] 
and between 2.2 and 2.6 cm for the synthetic fuel-rich furnace experiments [22] 
using the expression [24]: 


5 = 


1.5L 


Re v2 Sc va (17) 

All terms are the same as in Equation 1 1 . Since the furnace tube inner radius is 
1.1 cm, boundary layers of the sizes calculated above using Equation 17 are 
physically impossible. For these cases, fluxes of volatile species were calculated 
by restricting the boundary layer to half the furnace tube radius, 0.5 cm, using 
the following expression [24]: 


5 


( 18 ) 



Again all terms are defined as in Equation 12. Restricted boundary layers result 
in measured fluxes higher than those calculated using Equation 12. 

The calculated fluxes were compared to experimentally measured fluxes. 
Chemical models for volatilization were chosen which best represented the 
experimental data. These models were then used to extrapolate from furnace 
test conditions to HPBR conditions using expressions such as Equation 16. 

RESULTS 

The calculated boundary layer controlled fluxes are listed in Table 4. 
Contributions of each species to the total flux are specified. Calculated fluxes 
are compared to experimental results under each test condition. The magnitude 
of the fluxes as well as the pressure and temperature dependence are used to 
evaluate whether the boundary layer calculations adequately describe the 
measured fluxes and which Si-O-H species are likely contributing to the Si0 2 
volatility/SiC recession. Calculated Langmuir fluxes are several orders of 
magnitude larger than measured fluxes in all cases, demonstrating that a 
gaseous boundary layer acts as a barrier to volatilization. 

Synthetic Fuel-Lean Conditions 

Experimental results [13] as well as calculated fluxes for synthetic fuel- 
lean furnace conditions are shown in Figure 3. Calculated boundary layer 
controlled fluxes of Si(OH )4 based on Hashimoto’s [7] measured thermodynamic 
data agree with experimentally measured fluxes. Data set K indicates that 
SiO(OH) 2 should also be forming in measurable quantities, whereas data set A 



predicts negligible vapor pressures of this volatile species. Use of data set K in 
the flux calculations results in a larger overall flux of volatile species as well as a 
larger temperature dependence than that found experimentally. Thus, the 
synthetic fuel-lean furnace experimental results agree best with boundary layer 
controlled fluxes of an Si(OH) 4 volatile species alone. 

Fuel-Lean Burner Ria 

Experimental results [1] and calculated boundary layer controlled fluxes 
for HPBR fuel-lean conditions are shown in Figure 4. Data set K predicts that 
SiO(OH) 2 will be important in addition to Si(OH) 4 while data set A predicts that 
only Si(OH) 4 is important under these conditions. The experimental data fall 
between the predictions in both magnitude and activation enthalpy. This 
suggests that SiO(OH) 2 is present, but in quantities intermediate to those 
predicted by Krikorian and Allendorf. 

The pressure exponent for the volatilization rate in the HPBR fuel-lean 
case was found to be 1.5010.13. While this is in exact agreement with the 
predicted pressure dependence in the case of Si(OH) 4 formation (Equation 16), 
this result is not consistent with the presences of appreciable amounts of 
SiO(OH) 2 . Therefore, the pressure dependence indicates that Si(OH) 4 is the 
predominant volatile species in fuel-lean combustion conditions. 

Synthetic Fuel-Rich Conditions 

Experimental results [22] and calculated fluxes for synthetic fuel-rich 
furnace conditions are shown in Figure 5. Calculated boundary layer controlled 
fluxes of SiO agree best with both the magnitude and temperature dependence 



of the experimental data. While Si(OH) 4 is predicted to be more important than 
SiO, and the magnitude of the calculated volatility rate based on both species fits 
the data quite well, the temperature dependence of the sum of these fluxes (220 
kJ/mol) is not large enough to adequately describe the data (473±195 kJ/mol). 
Data set K predicts that the formation of SiO is negligible but that SiO(OH) and 
SiO(OH) 2 dominate the volatility of Si0 2 under these conditions. In contrast, 
data set A predicts that Si(OH) 4 is the primary volatile species at 1500K, SiO is 
the primary species at 1700K, while SiO(OH) 2 is only found in small amounts at 
1700K. In summary, the synthetic fuel-rich furnace experimental results agree 
best with boundary layer controlled fluxes of an SiO volatile species alone, 
although silicon hydroxide and oxyhydroxide species may form in measurable 
amounts. 

Fuel-Rich Burner Rio 

Finally, experimental results [1] and calculated boundary layer controlled 
fluxes for HPBR fuel-rich conditions are shown in Figure 6. Data set K predicts 
that Si(OH) 4 , SiO(OH) 2 , and SiO(OH) are all important. Si(OH) 4 is most 
important at 1500K while SiO(OH) is most important at 1700K. On the other 
hand, data set A predicts that only Si(OH) - 4 and SiO are important under these 
conditions. SiO is predicted to become significant only at 1700K. The 
experimental data agree quite well in magnitude with the predictions of data set 
K although the temperature dependence of the measured values is significantly 
less than the predicted values. The magnitude of the measured flux is 
significantly higher than that predicted using data set A which predicts the 



oxyhydroxides are unimportant. This suggests that one or both oxyhydroxide 
species are important under these conditions. 

The pressure exponent for the volatilization rate in the HPBR fuel-rich 
case was found to be 1.74±0.20. Because multiple volatile species are 
expected to be present in this environment, it was attempted to isolate the 
pressure effects of Si(OH) 4 from the remaining species. This is shown in Figure 
7a. The flux due to Si(OH)4 at 165 IK and several pressures was calculated 
using Hashimoto’s data [7]. This calculated value was subtracted from the 
measured value. The remaining pressure dependence of 2.03 was attributed to 
other volatile species. Several observations can be made from this plot. First, 
Si(OH) 4 makes up only a small fraction of the volatile products formed under 
these conditions And second, the pressure exponent of 2.03 necessarily implies 
that the unknown volatile species have a stronger dependence on water vapor 
pressure than Si(OH) 4 . This contradicts the previous conclusion that the 
unknown volatile species are SiO(OH) and/or SiO(OH) 2 which have flux related 
pressure exponents of -1/4 and 1/2 respectively. The expected pressure 
exponent for Si2(OH)6 flux from Equation 10 is 2.0. Note also that this molecule 
would likely be formed in fuel-rich environments and not in fuel-lean 
environments due to the oxygen generated as a product of Reaction 1 0. While 
Krikorian [11] has calculated free energy functions for this molecule, no enthalpy 
data are available. Estimation of the vapor pressure of this species may be 
possible, but is beyond the scope of this paper. 



The reaction enthalpy for Si0 2 volatility in the HPBR fuel-rich case was 
found to be 159+9 kJ/mol. Again, because multiple volatile species are expected 
to be present in this environment, it was attempted to isolate the temperature 
effects of Si(OH) 4 and SiO from the remaining species. This is shown in Figure 
7b. The flux due to Si(OH) 4 and SiO at 1500, 1600, and 1700K was calculated. 
This calculated value was subtracted from the measured value. The remaining 
flux hqd a reaction enthalpy of 170 kJ/mol which can be attributed to the 
formation of other volatile species from Si0 2 . It can be seen from Table 1b that 
the enthalpies of reaction for both SiO(OH) 2 and SiO(OH) formation are much 
higher than 170 kJ/mol. This again indicates that some other volatile species is 
needed to explain the high pressure fuel-rich results. While the reaction 
enthalpy of Si 2 (OH)6 by reaction 10 is unknown, entropy considerations 
(Reaction 10 results in a net loss of vapor molecules) necessitate that the 
reaction enthalpy be relatively low compared to reactions 6 and 9 [25]. This is in 
qualitative agreement with the fuel-rich HPBR results. 

Extrapolation of Laboratory Furnace Results to HPBR Conditions 

The simple volatilization models developed for the synthetic combustion 
gas furnace tests and the experimental variables of pressure and velocity are 
now used to extrapolate to the volatilization rates under HPBR conditions. 
Comparisons to measured HPBR volatilization rates are made. These results are 
shown in Figures 8a and 8b. 

For the fuel-lean case assuming Si(OH) 4 formation only, Equation 1 6 was 
used to extrapolate volatilization rates from 50% H 2 0 at 1 atm total pressure and 
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a gas velocity of 4.4 cm/sec to 12.3% H 2 0 (for .J>=0.9) at 6.3 atm total pressure 
and a gas velocity of 20 m/sec. This extrapolation is labeled (1) in Figure 8a. 
The extrapolation somewhat over-predicts the magnitude of the volatility and 
under-predicts the temperature dependence. Nevertheless, the extrapolation 
gives a fair prediction of the order of magnitude of the volatility rates measured in 
the HPBR in fuel-lean conditions. 

Previous discussion has shown that the synthetic fuel-rich furnace tests 
are best modeled by SiO volatility, but this model is inadequate to explain the 
HPBR results. The following discussion shows the shortfall in predicted volatility 
when using the one atmosphere model to predict volatility at the conditions of the 
HPBR. In the fuel-rich environment, the reducing species present in the / 
combustion environment, CO and H 2 , play a significant role in the silica volatility. 
The dependence of SiO on the partial pressures of combustion gases similar to 
Equation 16 must be formulated. The reactions for SiO formation given by 
Equations 4 and 5 have therefore been summed for this simple model of the 
combustion environment: 

2Si0 2 + CO(g) + H 2 (g) = 2SiO(g) + 0O 2 (g) + H 2 0(g) (19) 

The flux of SiO(g) for this model can then be expressed as: 

( 20 ) 

and since the partial pressures of the gaseous; components at a constant <j> 
scale with each other as total pressure changes: 
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Jsio 06 ~\i 2 ~ at constant <> (21) 

P total 

At constant <|> and gas velocity, the SiO(g) flux should thus decrease as the total 
pressure is increased. This behavior differs from that of Si(OH) 4 and SiO(OH) 2 . 

Equation 20 was used to extrapolate volatilization rates based on SiO 
volatility from a synthetic 4> of 1.5 at 1 atm total pressure and a gas velocity of 0.4 
cm/sec to <t>=1 .8 at 6.3 atm total pressure and a gas velocity of 20m/sec. This 
extrapolation is labeled (2) in Figure 8b. While the formation of SiO(g) explains 
the results for the furnace testing in synthetic fuel-rich conditions, it is clearly 
inadequate to explain the volatility observed in fuel-rich conditions at the higher 
pressures in the HPBR. 

Since the thermochemical data predict that Si(OH) 4 is also present in fuel- 
rich environments, the synthetic fuel-lean volatility measured at 1 atm was used 
to extrapolate Si(OH )4 formation to the HPBR fuel-rich conditions. Equation 16 
was used to extrapolate volatilization rates from 50% H 2 0 at 1 atm total pressure 
and a gas velocity of 4.4 cm/sec to 10.3% H 2 0 (for <J>=1.8) at 6.3 atm total 
pressure and a gas velocity of 20 m/sec. This extrapolation is labeled (3) in 
Figure 8b. The extrapolated values for SiO and Si(OH) 4 , curves (2) and (3), are 
then summed and shown as curve (4) in Figure 8b. This extrapolation of SiO 
and Si(OH) 4 from furnace testing to the HPBR fuel-rich conditions under-predicts 
the magnitude of the volatility of silica under these conditions. This extrapolation 
again shows that other volatile silicon-containing species, such as SiO(OH) 2 , 
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SiO(OH), and/or Si 2 (OH) 6 may be important in fuel-rich combustion environments 
at high pressures. 

DISCUSSION 

A gaseous diffusion model has been used to describe the volatility of Si02 
scales formed on SiC under combustion conditions. The temperature, pressure, 
and velocity dependencies as well as the overall magnitude of the volatilization 
rate conform to this basic model. However, exact agreement between the 
calculated and measured fluxes is not always obtained. Concerns with the 
deficiencies of the input thermodynamic data and the gaseous diffusion model 
itself are now discussed. 

Assessment of Thermodynamic Data 

The accuracy of prediction based on the formation of a particular volatile 
species depends on the accuracy of the thermochemical data. Much of the 
thermochemical data for the Si-O-H system has not been measured 
experimentally or is incomplete. In this section, we attempt to evaluate the 
accuracy of the experimental data with reference to our own laboratory 
observations. Then, in turn, we attempt to evaluate the estimated and calculated 
data with reference to any available experimental data. Clearly, additional 
experimental measurements of thermochemical data for the silicon oxyhydroxide 
species would be helpful to characterize the volatilization/recession processes of 
SiCVSiC system in combustion environments. 
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SifOH)* Comparison of furnace tests to Hashimoto’s [7] experimentally- 
determined thermochemical data 

The vapor pressures of Si(OH) 4 predicted from Hashimoto’s transpiration 
studies are consistently higher than those observed in our laboratory (by less 
than an order of magnitude). Our results are probably lower due to partial 
suppression of Si 02 volatility from the fused quartz reaction tubes used in our 
studies. This mechanism of partially suppressed volatility should not affect the 
HPBR results in our test, where no other Si-base materials were present. 

Si(OH) 4 : Comparison of estimated and calculated thermochemical data to 

measured thermochemical data. 

The free energy of formation for Si(OH) 4 at 1600K from Hashimoto’s [7] 
experimental study as well as Krikorian’s [11] and Allendorf’s [20] calculated 
values are shown in Table 5. All three values are in fairly good agreement with 
each other. Krikorian's values slightly under-predict the stability of Si(OH) 4 
whereas Allendorf’s values slightly over-predict Si(OH) 4 stability. Hashimoto’s 
experimentally determined values were used in all models in this study. 

SiO(OH )2 and SiO(OH): Comparison of estimated and calculated 

thermochemical data to experimental observations. 

Qualitative observations in our laboratory [8] using a specialized molecular 
beam mass spectrometer for the reaction of silica and water vapor at 1 300°C 
and 1 atmosphere total pressure support the conclusion that both Si(OH) 4 and 
SiO(OH) 2 should be observed in appreciable quantities when SiC is present in 
combustion environments. The intensity of SiO(OH) 2 was observed in quantities 
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ten times less than Si(OH) 4 under conditions where Krikorian’s data [11] predict 
a vapor pressures of SiO(OH) 2 three times that of Si(OH) 4 and Allendorf’s data 
[20] predicts SiO(OH) 2 vapor pressures about thirty times less than those of 
Si(OH) 4 . Hashimoto [7] observed a square dependence of volatile species on the 
water vapor partial pressure for temperatures between 1100 and 1500°C under 
similar conditions to those studied using mass spectrometry, indicating that 
Si(OH) 4 is the primary volatile species. The combination of these results 
indicates that SiO(OH) 2 is present in silica-containing combustion environments 
in measurable quantities but in amounts about midway between the predictions 
of Krikorian’s and Allendorf’s data. At temperatures higher than 1300°C, 
SiO(OH) 2 will become more significant relative to Si(OH) 4 due to the difference in 
the temperature dependence of the formation reactions. At higher pressures 
SiO(OH) 2 will become less important relative to Si(OH) 4 due to its smaller 
pressure dependence. 

SiO(OH)2 and SiO(OH): Comparison of estimated and calculated 

thermochemical data to measured thermochemical data. 

The free energy of formation of SiO(OH) 2 and SiO(OH) from Krikorian’s 
[11] and Allendorf’s [20] calculated values as well as the experimentally 
determined values of Hildenbrand and Lau [9,10] at 2065K are shown in Table 6. 
This temperature is the only temperature at which experimentally determined 
data are available. Hildenbrand and Lau’s values were obtained using Knudsen 
effusion mass spectrometry. Krikorian’s estimates are quite close to the 
measured values for both molecules, but slightly over-predict the stability of the 
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oxyhydroxides. Allendorf’s calculated free energies significantly under-predict the 
stability of both oxyhydroxide molecules. There is some lingering uncertainty 
regarding the experimentally determined values for SiO(OH) due to the complex 
fragmentation patterns occurring during mass spectrometry [20,10]. Allendorf 
[20] and Darling&Schlegel [21] calculated stability of SiO(OH) are in better 
agreement with Hildenbrand and Lau’s revised values for SiO(OH) [10] than 
Krikorian’s estimate. This revision results in a higher temperature dependence 
for reaction 9. 

Flow Conditions in the HPBR 

While the gas flow conditions in the furnace tests are well defined and fall 
well within the laminar flow regime, the gas flow conditions in the HPBR are more 
complicated and require some discussion. First, Reynold’s numbers (Re D ) 
calculated for nitrogen flow in the burner rig are about 8x1 0 4 , indicating turbulent 
flow conditions [18]. Re x calculated for a flat plate geometry, similar to a SiC 
coupon in the burner rig, is on the order of 1 xIO 4 . For Re x less than 10 5 , flat 
plate samples within the turbulent flow of the burner rig should develop a 
laminar flow boundary layer given a well defined velocity profile [18]. However, 
since the burner rig essentially has the geometry of a short pipe (length to 
diameter ratio is less than 10), entry effects are important and a well defined 
velocity profile is not yet established as the flow impinges on the SiC sample 
under test [14]. Use of either laminar or turbulent flow equations under these 
conditions may not apply. Finally, when turbulent flow is established, the Re 
exponent in the flux equation is 0.8 [17]. Thus, the velocity exponent is expected 
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to vary between 0.5 for laminar flow and 0.8 for turbulent flow. The regression of 
the SiC recession rates in the HPBR yielded velocity exponents of 0.50±0.16 for 
the fuel-lean case and 0.69±0.32 in the fuel rich case. The fuel-lean velocity 
exponent is in excellent agreement with the predicted laminar flow velocity 
exponent of 0.5. The fuel-rich velocity exponent, while not statistically 
significantly different from 0.5, may indicate that laminar flow equations do not 
apply under these conditions. This may explain the difficulty encountered in 
modeling the fuel-rich volatility chemistry using laminar flow equations. 

Volatile Sp ecies Found for SiC in Fuel-Lean Combustion Conditions 

The volatilization and recession process for the Si0 2 /SiC system in fuel- 
lean combustion conditions has been well characterized in this study. A 
summary of the identification of the volatile species based on magnitude of flux 
as well as pressure and temperature dependence is found in Table 7. Laminar 
flow analyses appear to adequately model the process. Volatilization of silica 
can be adequately explained by the formation of Si(OH) 4 alone. There is some 
evidence that SiO(OH) 2 is also formed under these conditions in measurable 
quantities. The importance of SiO(OH) 2 will increase with temperature and 
decrease with pressure. 

Volatile Sp ecies Found for SiC in Fuel-Rich Comhustion Conditions 

The volatilization and recession process for the Si0 2 /SiC system in fuel- 
rich combustion conditions is still not completely understood. The likely identity 
of the volatile species formed under these conditions based on the magnitude of 
the flux, the temperature dependence and the pressure dependence is also 
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summarized in Table 7. Volatilization at low pressures in the fuel-rich 
combustion environment (1 atm) can be attributed to the formation of SiO. Since 
the flux of this species decreases as total pressure is increased, it is not 
surprising that completely different volatile species are needed to explain the 
results at higher pressures. It is therefore difficult to model high pressure fuel- 
rich combustion environments in a one atmosphere furnace test. In addition, 
laminar flow may not be operative under the high pressure conditions studied 
here. The high pressure fuel-rich volatility mechanism is best described, given 
the current thermochemical data, by the formation of Si(OH) 4 , SiO(OH), and 
SiO(OH) 2 . SiO(OH) is expected to increase with importance as the temperature 
is raised and as the combustion environment becomes more fuel-rich. The 
Si 2 (OH) 6 species satisfies the pressure dependence of the Si 02 volatilization rate 
observed here, but little is known about this molecule. Despite the uncertainty in 
the chemical model for silica volatility in fuel-rich conditions, the multiple linear 
regression model (volatility rate expressed as a function of pressure, 
temperature, and velocity) obtained in Part I of this paper is useful to estimate 
the volatility rate in fuel-rich combustion environments. 

Conclusions 

Water vapor is responsible for the degradation of SiC in the fuel-lean 
environment. Silica volatility in fuel-lean combustion environments is attributed 
primarily to the formation of Si(OH) 4 (g) with a small contribution of SiO(OH) 2 (g). 
Water vapor in combination with the more reducing environment contributes to 
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the higher rates of degradation of Si02 in fuel-rich combustion conditions. Silica 
volatility in a fuel-rich combustion environment is best described by the formation 
of SiO(g) at one atmosphere total pressure and greater contributions from the 
formation of Si(OH) 4 (g), SiO(OH) 2 (g), and SiO(OH)(g) at higher pressures. The 
molecule, Si 2 (OH)e(g), would provide the correct pressure dependence of the 
volatility of Si0 2 under fuel-rich conditions, however, complete thermochemical 
data are unavailable for this species. Clearly the need exists for more 
experimental measurements of thermochemical data for the Si-O-H system to 
aid in the modeling of silica volatility in combustion environments. 
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Table la 

Pressure Dependence of Si-O-H(g) Species 


Reaction number 

Species 



6 




7 

E25239H 

n=2 

m=1.5 

8 

Si 2 0(0H) 6 

n=3 

m=2.5 | 

9 

SiO(OH) 

n=0.25 

m=-0.25 1 

10 

rSi 2 (OH) 6 

n=2.5 

m=2 | 


Table 1b 

Temperature Dependence of Si-O-H(g) Species, 
Reaction Enthalpy for Si0 2 + H 2 0 = Si-O-H(g) 


Species 

Source for 

Thermochemical Data 

Reaction 
Enthalpy, kj/mol 

Si(OH) 4 

Hashimoto* 

57 


krikorian 

32 


Allendorf 

60 | 

SiO(OH) 2 

Hildenbrand* 

— 260 


Krikorian 

247 


Allendorf 

~ 346 

SiO(OH) 

[Hildenbrand* 

530 


Hildenbrand revised* 

670 ' 


Krikorian 

524 | 


Allendorf 

_ — 


Darling & Schleael 

721 


experimentally determined, all others estimated or calculated 
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Table 2 

Parameters used for calculation of boundary layer controlled fluxes. 
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Table 3 

Thermochemical Data Sets Used for Flux Calculations 







Table 4 

Calculated Boundary Layer Controlled Fluxes 
(mg/cm 2 h) 





32 

Table 5 

The free energy of formation for Si(OH) 4 
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Table 6 

The free energy of formation for SiO(OH) and SiO(OH) 2 


AG°f (2065K) kJ/mol 


Hildenbrand & Lau 

Krikorian 

Allendorf 


KzzmBLmm 

-450.78 

-254.39 

1 35SJJI 

-525.66 

-546.65 

-454.38 


‘revised value 
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Table 7 

Probable Identification of Volatile Species in Combustion Environments 



Criteria for choice of 
volatile species 

Most likely volatile species 

Synthetic fuel-lean 

magnitude of flux 
T dependence of flux 

Si(OH) 4 " 

Si(OH) 4 

Synthetic fuel-rich 

magnitude of flux 
T dependence of flux 

SiO 

SiO 

HPBR fuel-lean 

magnitude of flux 
T dependence of flux 
P dependence of flux 

Si(OH) 4 , SiO(OH) 2 
Si(OH) 4 , SiO(OH) 2 
Si(OH) 4 

HPBR fuel-rich 

magnitude of flux 
T dependence of flux 
P dependence of flux 

SiO(OH), SiO(OH) 2 , Si(OH) 4 
Si 2 (OH) 6 (?), Si(OH) 4 
Si 2 (OHk Si(OH) 4 



in 

co 


0 

o> 

CD 

CO 

E 

O 

-C 

O 


•O 

0 

0 

3 

0 

0 

’o 

0 

Q. 

CO 

0 

13 

O 

> 


0 

c 

o 

'55 

0 

£ 

Q. 

X 

Lii 

15 

o 

E 

0 

-C 


•S o 


“O 

c 

o _ 

Q_ a) 
Q-JC 

< I- 


0 

0 


o "E 

*o 
__ c 
=3 iS 
0 cn 
£ 0 

Q_ -C 


O 

0 
i— , 

3 

CD 


0 
Q_ 

E 

.2 x - ® t£ 

© 2 J ^ 

«0 9_ 

o CD f- 


irJ <D . 
c O c/5 LU 

o w JS T 

J O Q^, 

3 >*© ir 

— q 

ID 5 fl , 

£ 85s 

— • * 

\- 

* 


f IJo 

^ w r + 

E E 

CD 


<D 

sz 

O 


© I- 


O - 
* 


<3 m 


+ 

< 


~ •£ co 
■a - cn ii 
© aj c\j _ 

S€ *■" 


M © m U 


© <15 
O «g .. 
z r E 


O) 

c 


o = 


c 
o 

0 — -C w 
0-0-0 
£ g 0 ° 

9* o JQ 

t- (0 

® =6 Vi 5 

p o © = 
jr o -c o 

c C w © 

© O ffl £ 

<D "ifl "O 

© 3 © 

4= -Q C 

© E © 

£ 8 ^ 

*o © ’” 

.1 § is 

iS w p 
Q. © CC 
co W 

O 2.52 

CO © © 

3 c 5 
-f £(D 

,E U Q) £ 

5 • DC H- 
O CO ^ . 

0 c 


0 


o 

0 

0 


O 

0 

c 

g 

V- 

0 

D 

cr 

0 

>s 

o> 

l. 

0 

c 

0 


o 0 

h- nr 


_ g © 

I -s J ! 

I- > LU .E 






































37 


* 


List of Figures 

Figurel. Equilibrium calculation of gas products from combustion of Jet A fuel 
and air as a function of equivalence ratio, ty. 

Figure 2. Paralinear oxidation kinetics for SiC in water vapor. Rate constants 
are typical of those found for 50% H 2 O/50% 0 2 at 1 atm total pressure, gas 
velocity of 4.4 cm/s, and temperature of 1200°C. A) Weight change kinetics. B) 
Dimensional change kinetics. 

Figure 3. Linear weight loss rates for SiC in synthetic fuel-lean conditions. The 

dashed line shows the measured rates while the solid lines show calculated 
rates. 

Figure 4. Linear weight loss rates for SiC in fuel-lean conditions in the HPBR. 

The dashed line shows the measured rates while the solid lines show calculated 
rates. 

Figure 5. Linear weight loss rates for SiC in synthetic fuel-rich conditions. The 

dashed line shows the measured rates while the solid lines show calculated 
rates. 

Figure 6. Linear weight loss rates for SiC in fuel-rich conditions in the HPBR. 

The dashed line shows the measured rates' while the solid lines show calculated 
rates. 

Figure 7. A) Pressure dependence (at about 1650K) and B) Temperature 

dependence of the flux of the unknown Si-O-H species from SiC in high-pressure 
fuel-rich conditions. 
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Figure 8. Comparison of weight loss rates for SiC measured in the HPBR to 
those extrapolated from furnace results based on A) Si(OH) 4 in fuel-lean 
conditions, and B) SiO + Si(OH) 4 in fuel-rich conditions. 
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